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Abstract 

The standing quasi modes in a cyhndrical incompressible flux tube with magnetic twist 
that undergoes a radial density structuring is considered in ideal magnetohydrodynamics 
(MHD). The radial structuring is assumed to be a linearly varying density profile. Using 
the relevant connection formulae, the dispersion relation for the MHD waves is derived and 
solved numerically to obtain both the frequencies and damping rates of the fundamental and 
first-overtone modes of both the kink (to = 1) and fluting (to = 2, 3) waves. It was found that 
a magnetic twist will increase the frequencies, damping rates and the ratio of the oscillation 
frequency to the damping rate of these modes. The period ratio P1/P2 of the fundamental 
and its first-overtone surface waves for kink (to = 1) and fluting (to = 2,3) modes is lower 
than 2 (the value for an untwisted loop) in the presence of twisted magnetic field. For 
the kink modes, particularly, the magnetic twists B^/B^ =0.0065 and 0.0255 can achieve 
deviations from 2 of the same order of magnitude as in the observations. Furthermore, for 
the fundamental kink body waves, the frequency bandwidth increases with increasing the 
magnetic twist. 
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1 Introduction 



Transverse oscillations of coronal loops were first identified by Aschwanden et al. (1999) and 
Nakariakov et al. (1999) using the observations of Transition Region and Coronal Explorer 
(TRACE). Nakariakov et al. (1999) reported the detection of spatial oscillations in five coronal 
loops with periods ranging from 258 to 320 s. The decay time was 14. 5 it 2. 7 min for an oscillation 
of 3.9ib0.13 mHz. Also Wang and Solanki (2004) described a loop oscillation observed on 17 April 
2002 by TRACE in 195 A. They interpreted the observed loop motion as a vertical oscillation, 
with a period of 3.9 min and a decay time of 11.9 min. Nakariakov et al. (1999) obtained that 
the amplitude of the oscillations decreases by more than 50 percent in several oscillation periods. 
To estimate the energy flux of the EUV kink oscillation observed by Nakariakov et al. (1999), we 
use the maximum kinetic energy flux of an oscillating loop given by ^{t: L) pv^^y./ {2'kRL)tb = 
jRpv'^g^^/TD. Where p, Umax, R, L and td are the mass density, the peak loop velocity, the loop 
radius, the loop length, and the damping time of oscillation, respectively. Taking R = 2 x 10^ 
km, p = 2 X 10~^^ gr cm^^ for a typical coronal loop in active region and using Umax = 47 km 
s~^ and Td = 14.5 min given by Nakariakov et al. (1999), we obtain the energy flux as 2 x 10^ 
erg cm~^ which is three orders of magnitude smaller than the radiated energy flux 10^ erg 
cm~^ s~^ in active regions (see Klimchuk 2006). Therefore, the energy flux in the EUV kink 
oscillations is unlikely to be sufficient to cover heating of coronal loops. 

Since the discovery of the coronal green line during the 1869 eclipse which identified as 
Fe XIV spectral line by Grotrian (1939), different theories of coronal heating have been put 
forward and debated (for reviews see, e.g., Hollweg 1991: Narain and Ulmschneider 1996; Walsh 
and Ireland 2003; Erdelyi 2004, 2005; Klimchuk 2006; Erdelyi and Ballai 2007; Taroyan and 
Erdelyi 2009). lonson (1978) was first to suggest that the resonant absorption of MHD waves in 
coronal plasmas could be a primary mechanism in coronal heating. Since then, much analytical 
and numerical work has been done on the subject. Rae and Roberts (1982) investigated both 
eikonal and differential equation approaches for the propagation of MHD waves in inhomogeneous 
plasmas. Hollweg (1987a,b) considered a dissipative layer of planar geometry to study the 
resonant absorption of coronal loops. Davila (1987) derived the resonant heating rate in the 
low-beta compressible resistive MHD approximation with the shear viscosity and found that the 
heating does not depend explicitly on the dissipation coefficients. Poedts, Goossens, and Kerner 
(1989, 1990) developed a finite element code to elaborate on the resonant absorption of Alfven 
waves in circular cylinders. 

Sakurai, Goossens, and Hollweg (1991a,b), Goossens, Hollweg, and Sakurai (1992), Steinolf- 
son and Davila (1993) did much work on resonant absorption. Ofman, Davila, and Steinolfson 
(1994), and Erdelyi and Goossens (1994) included viscous and resistive dissipations in their 
analysis and concluded that the heating rate due to shear viscosity is comparable in magnitude 
to the resistive resonant heating. Also, they concluded that the heating caused by compressive 
viscosity is negligible. Goossens, Ruderman, and Hollweg (1995) extended the analysis of Saku- 
rai et al. (1991b) on resonant Alfven waves in nonuniform magnetic fiux tubes for a static MHD 
equilibrium. They showed that the conservation law found in ideal MHD is valid in dissipative 
MHD. They derived analytical solutions to the dissipative MHD equations for the Lagrangian 
displacement and for the total pressure perturbation. Erdelyi and Goossens (1995) studied the 
heating of solar coronal loops by resonant absorption of Alfven waves in visco-resistive MHD. 
They pointed out that under solar coronal conditions the two dissipative mechanisms are both 
operational and that electrical resistivity appears to be slightly more important than viscosity. 
Erdelyi and Goossens (1996) studied the effect of an equilibrium flow on resonant absorption of 
linear MHD waves in compressible viscous cylindrical magnetic flux tubes. They showed that 
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the presence of the equiUbrium flow is very determinant for resonant absorption and signiflcantly 
affect the resonance absorption rate suffered by the incoming driving waves. Goossens, Andries, 
and Aschwanden (2002) used the TRACE data of Ofman and Aschwanden (2002) to infer the 
width of the inhomogeneous layer for 11 coronal loops. Rudcrman and Roberts (2002) solved 
the initial value problem for a zero-beta loop driven by a kink perturbation. They pointed out 
that the observed damping rate of the coronal oscillations is due to the resonant damping of the 
quasi-mode kink oscillations. They estimated the inhomogeneity length scale (thickness of the 
boundary layer) of the loop with the data of Nakariakov et al. (1999). Note that Nakariakov et 
al. (1999) showed that the damping rate of the global kink oscillation observed by the TRACE 
can be justified by the energy absorption due to anomalous viscosity. Whereas Ruder man and 
Roberts (2002) cleared that for large values of the viscous Reynolds number, the damping rate 
is independent of the viscous Reynolds number and proportional to the ratio of the thickness of 
the inhomogeneous layer and the radius of the tube. Van Doorsselaere et al. (2004a) used the 
LEDA code (Large-scale Eigenvalue solver for the Dissipative Alfven spectrum) to study the 
resistive absorption of the kink modes of cylindrical models. They concluded that, when the 
width of the nonuniform layer was increased, their numerical results differed by as much as 25 
per cent from those obtained with the analytical approximation. 

Safari et al. (2006) studied the resonant absorption of MHD waves in magnetized flux 
tubes with a radial density inhomogeneity. Using the approximation that resistive and viscous 
processes are operative in thin layers surrounding the singularities of the MHD equations, they 
concluded that as the longitudinal wave number increases, the maximum amplitude of the body 
eigenmodes shifts away from the resonant layer and causes a decrease in damping rates. 

An additional features of the flux tube is that of twist. There is observational evidence for 
the existence of twist in the solar atmosphere. The TRACE 171 A observations of a postflare 
loop system on 1999 conflrmed the existence of the twisted loops in the corona (see Aschwanden 
2005). Rotational movement along a loop observed by Chae et al. (2000) indicates the existence 
of twist. Aschwanden (2005) discusses observed noncoplanarity of loops which implies kinked 
field lines which, in turn, are also intimately related to twist. Besides, it is unfeasible that 
every fiux tube within the solar atmosphere is entirely twist free in spite of the random contin- 
uous motions observed at the footpoints (see Erdelyi and Carter 2006). Following Aschwanden 
(2005), the magnetic twist parameter B(f,/Bz can be estimated as B^/Bz = tanO. Here B^ and 
Bz are the azimuthal and axial components of the magnetic field, respectively. Also 9 is the 
shear angle between the untwisted and the twisted field line. The geometric shear angle 9, can 
observationally be measured in twisted coronal loops. 

Bennett, Roberts, and Narain (1999) examined the infiuence of magnetic twist on the modes 
of oscillations of a magnetic flux tube. They found that twist introduces an infinite band of body 
modes. Erdelyi and Fedun (2006) studied the wave propagation in a twisted cylindrical magnetic 
flux tube embedded in an incompressible but also magnetically twisted plasma. They found that 
increasing the external magnetic twist from to 0.3 caused an increase in the normalized periods 
of sausage MHD waves approximately by 1 — 2%. Erdelyi and Fedun (2007) extended the work 
of Erdelyi and Fedun (2006) to the case of compressible plasma. They found that increasing the 
internal twist may result in about up to 3 — 5% changes in periods for sausage modes. 

Erdelyi and Carter (2006) studied the propagation of MHD waves in a fully magnetically 
twisted flux tube consisting of a core, annulus and external region. They investigated their 
analysis by considering magnetic twist just in the annulus, the internal and external regions 
having straight magnetic fleld. Two modes of oscillations occur in this conflguration; surface 
and hybrid modes. They found that when the magnetic twist is increased the hybrid modes 
cover a wide range of phase speeds, centered around the annulus, longitudinal Alfven speed for 
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the sausage modes. 

Carter and Erdelyi (2007) investigated the oscillations of a magnetic flux tube configuration 

consisting of a core, annulus and external region each with straight distinct magnetic field in 
an incompressible medium. They found that there are two surface modes arising for both the 
sausage and kink modes for the annulus-core model where the monolithic tube has solely one 
surface mode for the incompressible case. Also they showed that the existence and width of an 
annulus layer has an effect on the phase speeds and periods. 

Carter and Erdelyi (2008) used the model introduced by Erdelyi and Carter (2006) to include 
the kink modes. They found for the set of kink body modes, the twist increases the phase 
speeds of the modes. Also they showed that there are two surface modes for the twisted shell 
configuration, one due to each surface, where one mode is trapped by the inner tube, the other 
by the annulus itself. 

Mikhalyaev and Solov'ev (2005) investigated the MHD waves in a double magnetic fiux tube 
embedded in a uniform external magnetic field. The tube consists of a dense hot cylindrical 
core surrounded by a co-axial shell. They found two slow and two fast magnetosonic modes can 
exist in the thin double tube. 

Verwichte et al. (2004), using the observations of TRACE, detected multimode oscillations 
for the first time. They found that two loops are oscillating in both the fundamental and the 
first-overtone standing kink modes. According to the theory of MHD waves, for uniform loops 
the ratio of the period of the fundamental to the period of the first overtone is exactly 2, but 
the ratios found by Verwichte et al. (2004) are 1.81±0.25 and 1.64zb0.23. However, these values 
were corrected with the improvement of the observational error bars to 1.82ib0.08 and 1.58ib0.06, 
respectively, by Van Doorsselaere, Nakariakov, and Verwichte (2007). Also Verth, Erdelyi, and 
Jess (2008) added some further corrections by considering the very important effects of loop 
expansion and estimated a period ratio of 1.54. All these values clearly differ from 2. This may 
be caused by different factors such as the effects of curvature (see e.g. Van Doorsselaere et al. 
2004b), leakage (see De Pontieu, Martens, and Hudson 2001), density stratification in the loops 
(see e.g. Andries et al. 2005; Erdelyi and Verth 2007; Karami and Asvar 2007; Safari, Nasiri, 
and Sobouti 2007; Karami, Nasiri, and Amiri 2009), magnetic field expansion (see Verth and 
Erdelyi 2008; Ruderman, Verth, and Erdelyi 2008; Verth et al. 2008) and magnetic twist (see 
e.g. Erdelyi and Fedun 2006, 2007; Erdelyi and Carter 2006; Karami and Barin 2009). 

Karami and Barin (2009) studied both the oscillations and damping of standing MHD sur- 
face and hybrid waves in coronal loops in presence of twisted magnetic field. They considered a 
straight cylindrical incompressible flux tube with magnetic twist just in the annulus and straight 
magnetic field in the internal and external regions. They showed that both the frequencies and 
damping rates of both the kink and fiuting modes increase when the twist parameter increases. 
They obtained that the period ratio P1/P2 of the fundamental and first-overtone for both the 
kink and fluting surface modes are lower than 2 (for untwisted loop) in presence of the twisted 
magnetic field. Ruderman (2007) studied nonaxisymmctric oscillations of thin twisted magnetic 
tubes in a zero-beta plasma by taking into account the longitudinal density stratification. Using 
the asymptotic analysis, he showed that the eigenmodes and eigenfrequencies of the kink and 
fiuting oscillations are described by a classical Sturm-Liouville problem for a second-order ordi- 
nary differential equation. He also concluded that the results concerning nonaxisymmctric waves 
in twisted magnetic tubes obtained by Bennett et al. (1999) for incompressible plasmas can be 
applied to global nonaxisymmctric waves, i.e. kink and fluting modes, in coronal loops despite 
that the coronal plasma is a low-beta plasma. This conclusion was also in a good agreement 
with Erdelyi and Fedun (2007). 

In the present work, our aim is to investigate the effect of twisted magnetic field on the 
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resonant absorption of standing MHD waves in the coronal loops to justify the rapid damping 
of oscillations and deviation of the period ratio P1/-P2 from 2 observed by TRACE. This paper 

is organized as follows. In Section 2 we derive the equations of motion, introduce the relevant 
connection formulae and obtain the dispersion relation. In Section 3 we give numerical results. 
Section 4 is devoted to conclusions. 

2 Equations of Motion 

The linearized MHD equations for an incompressible plasma are 

d5v ^ _V5p ^ 1 X 5B) X B + (V X B) X 5^} + -V^^v, (1) 
dt p Airp p 

f^ = Vx(ivxB) + — V^iB, (2) 

V • (5v = 0, (3) 

where Sv, (^B and 6p are the Eulerian perturbations in the velocity, magnetic field and ther- 
mal pressure, respectively; p, a, v and c are the mass density, the electrical conductivity, the 
viscosity and the speed of light, respectively. In the momcntinn Eq. (1) like Rudcrman and 
Roberts (2002) we write the viscous force in a simplified form vV^S^ instead of the classical Bra- 
ginskii's expression for the viscosity tensor in a magnetized plasma (Braginskii 1965). Because 
according to Ruderman and Roberts (2002), under typical coronal conditions, the coefficient of 
the shear viscosity is at least 10 orders of magnitude smaller than that of the compressional 
viscosity. However, in the problem of oscillations of coronal loops, dissipation is only impor- 
tant in an Alfvenic dissipative layer embracing an ideal resonant magnetic surface. Numerical 
studies by Ofman et al. (1994), and Erdelyi and Goossens (1994, 1995) have shown that in 
Alfvenic dissipative layers only the shear viscosity is significant, all other terms in Braginskii's 
tensorial expression being neglected. Equation (3) satisfies the incompressibility condition. The 
assumption of incompressibility reduces the direct full applicability to coronal loops except for 
kink modes observed by the TRACE as they are highly incompressible perturbations (see Carter 
and Erdelyi 2007). Recently Goossens et al. (2009) showed that in the thin tube approxima- 
tion neglecting contributions proportional to (/c^i?)^ then the frequency of the kink wave and 
its damping due to resonant absorption are the same in the three cases including a compress- 
ible pressureless plasma, an incompressible plasma and a compressible plasma which allows for 
MHD radiation. This is expected as kink modes are linearly incompressible. Note that here the 
energy equation is absent. Because for incompressible plasmas when sound speed (or adiabatic 
gas index) goes to infinity then the energy equation is decoupled from the system of linearized 
MHD equations. 

The simplifying assumptions are as follows. 

• The background magnetic field is assumed to be 

B = (0,^r,5^), 

where A, Bz are constant and the magnetic field is uniformly twisted inside and outside 
the tube (see Fig. 1). This is unphysical as r ^ 00, and caution has to be exercised. 



5 



Following Ruderman (2007), to satisfy the Shafranov-Kruskal stability criterion we assume 
that the azimuthal component of magnetic field is smaller than the axial component. This 
assumption is also compatible with observed weakly twisted coronal loops (see Erdelyi 
and Fedun 2006). The stability of a twisted magnetic flux tube has been studied in some 
details by Bennett et al. (1999), Ruderman (2007), and Carter and Erdelyi (2008). Note 
that the choice of B^j, = Ar implies that there is a constant longitudinal current density 
along the flux tube. 

Prom Safari et al. (2006) and Karami et al. (2009), the density profile is assumed to be 

f PU {r < Ri), 

Pir-) = l [^]iR-r) + Pe, iRi<r<R), 
[ Pe, {r > R), 

where R denotes the radius of the tube and Ri < R is the radius at which the resonant 
absorption occurs. Also pi and pe are the interior and exterior constant densities of the 
tube, respectively (see Fig. 1). Note that for resonant absorption a number of density 
profiles have been considered. For instance, Ofman et al. (1994, 1995) have taken p{r) = 
Pr + {1 — pr)e~^ with pr = Pe/ Pi- Here since we use the thin boundary approximation, we 
assume that the inhomogeneous layer and the resonance layer coincide (see Goossens et al. 
2009). Therefore any well behavior function of the density like that selected by Ofman et 
al. (1994, 1995) in the resonance layer, i.e. the region Ri < r < R where the singularity 
occurs, can be approximated with a linearly varying profile. 



The equilibrium condition, i.e. radial force balance equation ^ p(r) + 
gives the equilibrium plasma pressure p{r) as 

P(r) =Pq- 



B\r) 
Stt 



where po is the plasma pressure at the center of the tube. Since we have no discontinuity 
in the background magnetic field, then the thermal pressure remains continuous across 
r = R (see Fig. 1). 

• Tube geometry is a circular with cylindrical coordinates, (r, (j), z). 

• There is no initial steady flow over the tube. 

• Viscous and resistive coefficients, v and a respectively, are constant. 

• t-, (p- and z- dependence for any of the components (5v and 5B is exp{i(m^ + kzZ — ujt)}. 
Here kz = lir/L, L is the length of the tube, and Z = (1, 2, • • •), m = (0, 1,2,- • •) are the 
longitudinal and azimuthal mode numbers, respectively. 

In the remainder of this section the following steps are taken: a) in r < Ri and r > R, 
dissipative terms are neglected. Solutions of Eqs. (1) and (2) are obtained as in Bennett et al. 
(1999), Erdelyi and Carter (2006), and Karami and Barin (2009); b) in i?i < r < R, within 
which the resonant layer resides, we use the thin boundary (TB) approximation in which we 
assume that the inhomogeneous layer and the dissipative layer coincide (see Goossens et al. 
2009). This enables us to avoid solving the non-ideal MHD equations in the inhomogeneous 
layer. The jumps across the resonant layer are found by the prescriptions of Sakurai et al. 
(1991a); c) substituting the solutions (a) in jump conditions (b) gives an analytical expression 
for a dispersion relation to be solved for the frequencies and the damping rates. 
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2.1 Interior and Exterior Solutions 



Following Bennett et al. (1999), Erdelyi and Carter (2006), and Karami and Barin (2009), in 
the absence of dissipations, taking time derivative of Eq. (2) and substituting for dS\/dt from 
Eq. (1), the resulting equation yields to Bessel's equation for the Eulerian perturbation in total 
pressure SPt as 



where 



and 



d^ ^ 1 d 2^ 
dr^ r dr ^ ' 



6Pt = 0, 



5Pt = Sp + 



B(5B 

47r 



(4) 



(5) 



(6) 



{mA + k^B^). 



(7) 



Note that Eq. (4) is valid only outside the inhomogeneous layer. Equation (4) is same as the 
result exactly derived by Bennett et al. (1999), and Erdelyi and Carter (2006). Also Eq. (4) 
in its general form for a compressible plasma was derived by Erdelyi and Fedun (2007). Note 
that subscripts i (which correspond to the internal region) are replaced by e corresponding to 
the external region. Solutions of Eq. (4) for the interior region {r < R\) are: 



^ I Imimr), mf > 



surface waves. 



Jm{nir), nf = —mf > 0, body waves, 



(8) 



where Jm and Im are Bessel and modified Bessel functions of the first kind, respectively. In the 
exterior region {r > R), the waves should be evanescent. Solutions are 



6Pt = K^inier), ml > 0, 



(9) 



where is the modified Bessel function of the second kind. One must note that from Eq. (6), 

solutions (9) arc satisfied only for the frequencies u'^ > Lij\ + 



Sect. 3 confirms that this condition is hold for the obtained frequencies 



Our numerical results in 



2.2 Connection Formulae, Dispersion Relation and Damping 

According to the connection formulae given by Sakurai et al. (1991a), the jump across the 
boundary (resonance layer) for = —Svr/ico and SPt is 



[&] = -i 



ITT- 



1 93 



Ca, 



(10) 



[5Pt 



-ITT- 



1 AB^F 
A|27r/952 



Ca, 



(11) 
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where 



Ca = QboPt 27i~^''' 



dr 
m 
m5 



dr 



9B = --k,Ar. (12) 

r 

Note that a = R — Ri'is the thickness of the inhomogeneous layer and i?i < r a < -R is the radius 
at which the singularity occurs. Davila (1987) showed that in the resonance absorption, however, 
the damping rate is independent of the dissipation coefficient values. But the resonance layer 
width scales as 5a = + Estimates for typical coronal values suggest that resonant 

layers have thicknesses from 0.3 km to 250 km (Davila 1987). According to Ofman et al. (1994), 
and Erdelyi and Goossens (1994, 1995) for resonant absorption the main contribution comes not 
from bulk viscosity but from shear viscosity. If we use the Reynolds TZ = / (^^) = 560 

and Lundquist S = = 10^ numbers given by Ofman et al. (1994), and taking 

L = 109 X 10^ km, R/L = 0.01, a'/R = 0.08, pc/Pi = 0.1, and interior Alfvcn velocity va, = 2000 
km s~^ for a typical coronal loop, then one can get 5a — 85 km which is very close to the thickness 
of the inhomogeneous layer a 2± 87 km. Therefore we can use the thin boundary approximation 
which assumes that the thickness of the resonance layer is the same as the inhomogeneous layer 
width (see Goossens et al. 2009). 

Substituting the solutions of Eqs. (8)-(9) in jump conditions gives the dispersion relation as 

dQ{G:) + di{u:) = Q, (13) 

where a; = a; — ia, a is damping rate and 

_ miRr^{miRi) mFAR / meRK'^jmeR) mFA \ 



]j[ C'a f QBR Xe FAB^ / meRK^{meR) mFA \^ 

\A\p{rA)B^rA)Im{miR,)\ + Ai 27r V K^{m,R) ^27rDjr ^ > 



with 



A 



'J ~ n2 



Df - 4(FA/47r)2 ' 

Di = ^-Pi^'^ (16) 

where (j) stands for (i) or (e). Also a prime on and indicates a derivative with respect 
to their appropriate arguments. The results for the body waves are the same as Eqs. (14)-(15), 
except that Im is replaced by Jm everywhere. 

Equations (lO)-(ll) show that when the twist is absent, i.e. A = 0, the total Eulerian 
pressure will be continuous and only jumps across the boundary which is in agreement with 
Sakurai et al. (1991a). 

For the surface modes, solving Eq. (13) yields one mode whose frequency exists between oja-^ 
and coAe- But for the body modes, there is an infinite band of body modes with the frequencies 
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centered around ■ The existence of an infinite set of body waves in the presence of twist has 
been aheady introduced by Bennett et al. (1999). 

Note that the distinction between surface and body waves is not so important when we study 
the kink waves in slender magnetic tubes. Recently Goossens et al. (2009) concluded in the 
thin tube approximation, the kink MHD waves do not care about propagating (body wave) or 
evanescent (surface wave) behavior in the internal part of the flux tube. Although under coronal 
condition they are formally the body waves, they have typical properties of the surface waves. 
For example, the radial displacement takes its maximum value at the tube boundary. This is 
why Ruderman and Roberts (2002) suggested to call these waves the global kink waves. Also 
Erdelyi and Carter (2006) showed that under coronal conditions, for longer wavelength values 
the sausage surface mode only exists for smaller amounts of twist in the incompressible plasma. 

The numerical solution of the dispersion relation (13) yields to frequencies conmi and damping 
rates anmh which are characterized by a trio of wavenumbers (n, m and I) that actually count 
the number of nodes or antinodes along r-, 0- and z-directions, respectively. 

3 Numerical Results 

As typical parameters for a coronal loop, we assume L = 109 x 10^ km, i?/L = 0.01, a/R = 0.08, 
Pe/pi = 0.1, Pi = 2 X 10~^^ gr cm~^, = 100 G. For such a loop one finds va^ = 2000 km s~^, 
uj^. := ^ ^ 0.02 rad s"!. 

The effects of twisted magnetic field on both the frequencies a; and damping rates a are 
calculated by numerical solution of the dispersion relation, i.e. Eq. (13). The results are 
displayed in Figs. 2 to 9. Figures 2 to 7 show the frequencies, damping rates and also the ratio 
of the oscillation frequency to the damping rate of the fundamental and first-overtone I = 1,2 
kink (m = 1) and fluting (m = 2,3) surface modes versus the twist parameter, B^/Bz := ^ 
for relative inhomogeneous layer width a/R = 0.08. According to Carter and Erdelyi (2008), 
the kink (m = 1) speed is independent of the sound speed. Hence, the kink modes are highly 
incompressible, and this makes the limit of incompressibility be of great interest, not just for 
wave studies in the deeper part of the solar atmosphere (e.g. where the plasma-beta is high) 
but can also be directly applicable from the lower solar atmosphere to the corona. For all other 
modes, {m ^ 1), compressibility may be necessary. 

Figures 2 to 7 reveal that: i) the frequencies, the damping rates and the ratio Lo/a increase 
when the twist parameter increases. The result of u is in good agreement with that obtained 
by Carter and Erdelyi (2008). Also the behavior of co and a versus the twist parameter are in 
good concord with that obtained by Karami and Barin (2009). But, the numerical values of 
our damping rates for m =1, 2 and 3 are six, seven and eight orders of magnitude, respectively, 
greater than those obtained by Karami and Barin (2009). This is due to existence of the resonant 
absorption which is absent in their work. Note that the dissipative processes in Karami and 
Barin (2009) depend on the Reynolds TZ and Lundquist S numbers. They obtained the damping 
rates with TZ = 560 and S = 10^ given by Ofman et al. (1994). But in the resonance absorption, 
the damping rate is independent of the dissipation coefficient values and only the resonance 
layer width depends on TZ and S (see Davila 1987). ii) For m = 1 with increasing B^/B^, the 
frequencies, damping rates and their ratio, for instance com, CKiii and win/am, increase c±8-66, 
~6-50 and ~2-ll percent, respectively, compared with an untwisted loop, iii) For a given m and 
Bffi/Bz, when the longitudinal mode number, I, increases, both the frequencies and damping 
rates increase. But the ratio oj/a: for m = 1, for the twists lower than 0.003 increases and 
then decreases; for m = 2, always decreases; but for m = 3, for the twists lower than 0.115 
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decreases and then increases, iv) For a given I and B^/B^, when the azimuthal mode number, 
m, increases, the frequencies and damping rates increase but the ratio oj/a decreases. 

Here in our calculations, the sausage modes (m = 0) are absent. Because following Edwin 
and Roberts (1983), and Roberts, Edwin, and Benz (1984), the sausage modes have a cut-off 
longitudinal wave number l^., and they are only expected in fat and dense loops. According 
to Aschwanden (2005) the longitudinal wave number cut-off Ic for sausage-mode oscillations is 
expressed as a requirement of the loop length-to- width ratio L/{2R) as a function of the density 
contrast Pc/Pi between the external and internal loop densities. For instance, for a typical active 
region, loops which have a density contrast of the order of Pe/Pi ~ 0.1 — 0.5 would be required 
to have length-to- width ratios of L/{2R) 1 — 2. 

The period ratio P1/P2 of the fundamental and first-overtone, I = 1,2 modes of both the 
kink (m = 1), and fluting (m = 2, 3) surface waves versus the twist parameter is plotted in Fig. 
8. Figure shows that for both kink and fluting modes, the period ratio P1/P2 decreases when 
the twist parameter increases. For instance, P1/P2 decreases from 2 (for untwisted loop) and 
approaches below 1.6, 1.2 and 1.1 for m = (1, 2 and 3), respectively, with increasing the twist 
parameter. Note that when the twist is zero, the diagrams of P1/P2 do not start exactly from 
2. This may be caused by the radial structuring pe / Pi- But for the selected thin tube with 
R/L = 0.01, this departure is very small, 0(10-3) for kink (m = 1) and 0(10"^) for fluting 
(m = 2,m = 3) modes, and does not show itself in the diagrams (see McEwan et al. 2006). 
Figure 8 clears that for kink modes (m = 1), the ratio P1/P2 is 1.821 when i?^/i?2=0.0065, 
and is 1.584 for B^/Bz=0.0212. These are in good agreement with the period ratios observed 
by Van Doorsselaere et al. (2007), 1.82ib0.08 and 1.58ib0.06, respectively, deduced from the 
observations of TRACE. Van Doorsselaere et al. (2007) considered only the effect of density 
stratification and neglected the effect of loop expansion. It has been shown by Verth and Erdelyi 
(2008) that failure to take account of the flux tube expansion causes the density scale height to 
be overestimated by a factor of approximately 2. Verth et al. (2008) considered both effects of 
the longitudinal density stratification and magnetic field expansion and estimated the correct 
coronal density scale height. They reported a period ratio of 1.54 deduced from the observations 
of TRACE for the kink modes. Here in our numerical calculations, if we take 50/5^=0.0255 
and Pe/pi = 0.05 then we can obtain the period ratio Pi/P2=l-539 for the kink modes which 
is in good agreement with that reported by Verth et al. (2008). This shows that the magnetic 
twist can achieve deviations from 2 of the same order of magnitude as in the model of Verth 
et al. (2008) containing the longitudinal density stratification and magnetic field expansion. 
Note that for the kink modes, the only period ratio obtained by Karami and Barin (2009) was 
1.882 with i?0/i?2=O.OO65. Therefore considering the magnetic twist in presence of the resonant 
absorption not only improved the period ratio estimated by Karami and Barin (2009) but also 
justified another period ratio, 1.54, reported by Verth et al. (2008). 

Figure 9 displays the frequency band width, Aoj, including infinite set of the fundamental 
kink (m = 1) body modes versus the twist parameter. Figure 9 presents that Aoj increases when 
the twist parameter increases. This is in good agreement with the result obtained by Carter 
and Erdelyi (2008), and Karami and Barin (2009). 

4 Conclusions 

Resonant absorption of standing MHD surface and body waves in coronal loops in the presence 
of the twisted magnetic field is studied. To do this, a typical coronal loop is considered as a 
straight cylindrical incompressible flux tube with magnetic twist that undergoes a radial density 
structuring. Under the thin tube approximation, the thickness of the inhomogeneous layer is 
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same as the resonance layer width. Hence, the radial structuring is assumed to be a linearly 
varying density profile in the inhomogeneous layer. Using the relevant connection formulae, the 
dispersion relation is obtained and solved numerically for obtaining both the frequencies and 
damping rates of the fundamental and first-overtone kink and fluting modes. Our numerical 
results show that 

i) the frequencies and damping rates as well as the ratio of the oscillation frequency to the 
damping rate of both the kink (m = 1) and fluting (m = 2,3) surface waves increase when the 
twist parameter increases. For m = 1 with increasing B^/B^-, the ratio win/am changes from 
39.3 to 43.5 which approximately one order of magnitude is greater than the ratio reported by 
Nakariakov et al. (1999), Wang and Solanki (2004), and Verwichte et al. (2004) deduced from 
the TRACE data; 

ii) the period ratio PxjP^^ for both the kink (m = 1) and fluting {m = 2,3) surface modes 
is lower than 2 (for an untwisted loop) in the presence of the twisted magnetic field. The result 
of P1/P2 = 1.821, 1.539 for kink modes is in agreement with the TRACE observations; 

iii) Frequency bandwidth of the fundamental kink (m = 1) body modes increase when the 
twist parameter increases. 

Note that the main differences between the present work and the work of Karami and Barin 
(2009) can be summarized as follows: i) in our model, the magnetic field was assumed to be 
uniformly twisted inside and outside the tube. Whereas they considered only a uniformly twisted 
magnetic annulus inside the loop, ii) They considered the weak damping due to viscous and 
resistive dissipations which cannot justify the rapid damping of the coronal loops observed by 
Nakariakov et al. (1999), and Wang and Solanki (2004). In the case one includes the resonant 
absorption, the rapid damping can be justified, iii) Our model could justify the period ratios of 
kink modes reported by Van Doorssclaere et al. (2007), i.e. P1/P2 =1-82, 1.58, and Verth et al. 
(2008), i.e. P1/P2 =1.54, whereas Karami and Barin (2009) only justified P1/P2 =1.82. 
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Figure 1: Dimensionless equilibrium quantities including density p, Alfven velocity va, azimuthal 
component of the uniformly twisted magnetic field B^, thermal pressure p, and total pressure 

Ptot = P + S^/(87r) against fractional radius x = r/R. Auxiliary parameters are Pc/ Pi = 0.1, 
^ = 0.01, f3 = g2^°g^^ = 0.01. The parameter f5 has been taken from Smith, Tsiklauri, and 
Ruderman (2007). ' 
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Kink Surface Modes (m=l, 1=1) 




* z 

Figure 2: Frequency of the fundamental kink (m = 1) surface mode, its damping rate and the 
ratio of the oscillation frequency to the damping rate versus the twist parameter, S^/iJ^. The 
loop parameters are: L = 109 x 10^ km, R/L = 0.01, a/R = 0.08, pe/pi = 0.1, pi = 2 x 10"^^ 
gr cm~^, Bz = 100 G. Both frequencies and damping rates are in units of the interior Alfven 
frequency, — 0.02 rad s~^. 
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0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 




Figure 3: Same as Fig. 2, for the first-overtone kink (m = 1) surface modes. 
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Fluting Surface Modes (m=3, 1=1) 
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Figure 6: Same as Fig. 2, for the fundamental fluting (m = 3) surface modes. 
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Figure 7: Same as Fig. 2, for the first-overtone fluting (m = 3) surface modes. 
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Figure 8: The period ratio P1/P2 of the fundamental and its first-overtone surface modes versus 
the twist parameter for kink (m = 1) and fluting (m = 2, 3) modes. Auxiliary parameters as in 
Fig. 2. 
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Figure 9: Frequency bandwidth of the fundamental kink (m=l) body modes versus the twist 
parameter for the relative inhomogeneous layer width a/R = 0.08. Auxiliary parameters as in 
Fig. 2. 
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